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Abstract 

For AT > 2 supergravities, BPS black hole solutions preserving four supersymmetrics 
can be superposed linearly, leading to well defined solutions containing an arbitrary number 
of such BPS black holes at arbitrary positions. Being stationary, these solutions can be 
understood via associated non-linear sigma models over pseudo-Ricmaniann spaces coupled 
to Euclidean gravity in three spatial dimensions. As the main result of this paper, we 
show that whenever this pscudo-Ricmanniann space is an irreducible symmetric space 
<S/Sj*, the most general solutions of this type can be entirely characterised and derived 
from the nilpotent orbits of the associated Lie algebra g. This technique also permits the 
explicit computation of non-supersymmetric extremal solutions which cannot be obtained 
by truncation to Af = 2 supergravity theories. For maximal supergravity, we not only 
recover the known BPS solutions depending on 32 independent harmonic functions, but 
in addition find a set of non-BPS solutions depending on 29 harmonic functions. While 
the BPS solutions can be understood within the appropriate Af = 2 truncation of Af = 8 
supergravity, the general non-BPS solutions require the whole field content of the theory. 
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1 Introduction 



Extremal static black holes with identical charges do not interact. As first pointed 
out by Papapetrou and Majumdar, for every spherically symmetric extremal black hole 
of Maxwell-Einstein theory, one can arrive at associated multi-black hole solutions by 
means of linear 'superposition', that is, by substituting any harmonic function of the 
inverse radius in the potentials defining the stationary metric and the vector field [1, 2]. 
Physically, the stability (or 'staticity') of these solutions is ensured by the 'no force' 
property, whereby the gravitational attraction between any two charged black holes is 
exactly balanced by their mutual electromagnetic repulsion. On the other hand, con- 
sidering Maxwell-Einstein theory as the bosonic sector of M = 2 supergravity, one can 
understand this property mathematically from the fact that the corresponding solutions 
are BPS, as they preserve one-half (that is, four) of the supersymmetry charges. For 
general matter coupled M = 2 supergravity theories, there are multi-black hole solutions 
of BPS type for any black hole preserving four identical supersymmetry charges. Such 
solutions have been identified in a general way in [3], in the framework of the so-called 
attractor mechanism [4, 5]. 

However, the BPS multi-black hole solutions are not the most general solutions of this 
type. In this publication we characterise the general solutions of Papapetrou-Majumdar 
type for all models for which the scalar fields coordinatise a symmetric space. More 
precisely, our discussion applies to all models containing gravity coupled to abelian vec- 
tor fields and scalar fields parametrising a symmetric Riemannian space 04/^4 whose 
isometry group &4 acts faithfully on the vector fields. For these models, the station- 
ary equations of motion reduce to the ones of a non-linear sigma model defined over 
a pseudo-Riemaniann symmetric space <&/$)* coupled to three-dimensional Euclidean 
gravity [6, 7, 8], with G5 being a simple Lie group, and fi* a non-compact real form 
of its maximal compact subgroup. The derivation of Papapetrou-Majumdar type so- 
lutions then reduces to certain algebraic conditions [9] (see eqs. (3.2) and (3.5) of this 
paper). Using the formalism developed in [8], and especially the explicit relation be- 
tween extremal solutions and nilpotent orbits of the three-dimensional duality group (3, 
we are able to solve these algebraic equations in full generality. A crucial feature here is 
that, unlike supergravity duality groups in dimensions D > 4, the duality groups <5 of 
the three-dimensional theories also incorporate the Ehlers symmetry SX(2,IR), and thus 
the gravitational degrees of freedom. The importance of Lie algebra nilpotent orbits in 
the physics of extremal black holes was first pointed out in [10]. The case of minimal 
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supergravity in five dimensions has been studied in detail in this framework in [11]. 

In contrast to previous studies of BPS solutions, the analysis of [8] is solely based 
on exploiting properties of the (J5-Noether charge ^ and the associated BPS parameter 
c, which are defined in eqs. (2.5) and (2.8) below. As shown in [8], the Noether charge 
satisfies a cubic or quintic characteristic equation, cf. (3.13) and (3.14), respectively. 
The latter determines the BPS parameter as a function of the physical charges of the 
four- dimensional theory, namely the complex gravitational charge w = m + in (with m 
the mass and n the Taub-NUT parameter) and the complex (dyonic) electromagnetic 
charges Za = qA + ipA, viz. 

c = c(w,Z A ) (1.1) 

Solutions of Papapetrou-Majumdar type are superpositions of extremal black holes 
whose Noether charges satisfy the condition c = 0. These solutions come in two cate- 
gories. Those of the first category are \ BPS solutions of an appropriate truncation to 
a theory that constitutes the bosonic sector of some N = 2 supergravity theory coupled 
to n v vector multiplets. The | BPS solutions of the latter theory define multi-black hole 
solutions depending on 2 + 2n v harmonic functions (l + 2n v if one requires the solution to 
be asymptotically Minkowski). They define the most general solutions of Papapetrou- 
Majumdar type whenever the associated group $)* is non-semi- simple, i.e. possesses a 
£7(1) factor. In this case, the BPS parameter can be written in the form 



(k 


2 _ 


Zi\ 


2 )(k 


2 _ 




2 ) 




\w 


2 



(1.2) 



where z± and z 2 generalise the central charges associated to the electromagnetic charges 
to theories which are not necessarily super symmetric. It then follows that the extremality 
condition c = requires at least one of the charges to saturate the BPS bound \z\ = \w\. 
Consequently, the solution is \ BPS within an appropriate Af = 2 truncation. 

When Sj* is semi-simple, on the other hand, the BPS parameter c cannot be expressed 
in the simple form (1.2), but is a more complicated non-rational function of the physical 
charges. In this case, the condition c = does not necessarily imply saturation of any 
BPS bound. This is borne out by Papapetrou-Majumdar type solutions in the second 
category. As shown in [8] this can happen in particular for maximal M = 8 supergravity, 
or more generally for theories that can be obtained by dimensional reduction of five- 
dimensional supergravities on a circle. As we will show, such theories with n vector 
fields in four dimensions admit multi-black hole solutions depending on 1 + n harmonic 
functions (n for asymptotically Minkowski solutions) that do not define | BPS solutions 
of any N = 2 truncation of the theory. 
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Because they rely essentially only on properties of the three-dimensional duality 
groups (9, our results in principle allow a general classification. In order to bring out 
the main points as clearly as possible, we here restrict attention mostly to the example 
of maximal J\f = 8 supergravity [12]. The M = 4 supergravities will be discussed in 
detail in a companion publication [13]. The extremal spherically symmetric black holes 
of maximal supergravity can be understood within the simplified context of the so-called 
STU J\f = 2 supergravity [14]. From the point of view of the latter, there are three types 
of generic extremal spherically symmetric black holes: the \ BPS ones, the non-BPS ones 
for which the central charge vanishes on the horizon {Z = 0), and the non-BPS ones for 
which the central charge does not vanish on the horizon [Z ^ 0) [16, 17, 18]. Both | BPS 
and non-BPS black holes for which the central charge vanishes on the horizon correspond 
to | BPS black holes within maximal supergravity. Within the STU model such solutions 
are related by exchanging the vector fields belonging to the gravity multiplet with those 
in the matter multiplets, respectively. By contrast, the non-BPS solutions with Z ^ 
are still non-BPS when embedded in maximal supergravity. As we will explain in this 
paper, such solutions are associated to theories that define compactification on a circle 
of theories defined in 4 + 1 dimensions. For instance, pure gravity in five dimensions 
dimensionally reduced on a circle does admit extremal solutions, and it obviously defines 
a consistent truncation of any theory originating from five dimensions [11]. One then has 
a corresponding f)*-orbit of extremal solutions, which for the STU model corresponds to 
the non-BPS solutions with non-vanishing central charge on the horizon. 

The plan of the paper is as follows. We first review the general formalism [6, 8] and 
explain how the Papapetrou-Majumdar solutions can be read from the nilpotent orbits. 
Then we explain how one can understand the two categories of Papapetrou-Majumdar 
solutions within this formalism. We display in detail the spectrum of such solutions in 
maximal supergravity. In a subsequent publication by one of the authors [13], the case 
of M = 4 supergravity coupled to n vector multiplets will be discussed in detail. In 
particular, it will include explicit non-BPS multi-black hole solutions belonging to the 
second category (the would be non-BPS solutions with Z ^ within the STU model) 
depending on 6 + n independent harmonic functions. 

2 Extremal solutions and nilpotent orbits 

Stationary solutions of Einstein theory of gravitation in four dimensions coupled to 
abelian vector fields and scalar fields parametrising a symmetric space (55 4 /f) 4 can be 
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effectively obtained by solving the equations of motion of an associated non-linear sigma 
model coupled to gravity in three spatial dimensions (i.e. with Euclidean signature). 
When the isometry group ©4 acts on the abelian vector fields in a faithful representation 
[4 to define a non-linear symmetry of the equations of motion, one thus arrives at a three- 
dimensional non-linear sigma which unifies all the physical bosonic degrees of freedom 
into a single pseudo-Riemanian coset space (S/fj*, with (5 a simple Lie group and Sj* a 
non-compact real form of its maximal compact subgroup [6, 7]. 

To derive the equations of motion for a coset representative V in <5/Sj*, one decom- 
poses the Maurer-Cartan form V~ l dV into its coset and its f)* components, 

V~ 1 dV = Q + P , QEQ^efj* , P = P„dx» e q e f)* (2.1) 

where u, ■ • • = 1, 2, 3. The equation of motion for the scalars reads 

d*VPV _1 = (2.2) 

where * is the Hodge star operator associated to the three-dimensional Riemannian metric 
g^u. The Einstein equation is 

= Tr P^P V (2.3) 

Because in three dimensions the Riemann tensor is entirely determined by the Ricci 
tensor, the three-dimensional metric g^ is flat if the right hand side of (2.3) vanishes. 

Black hole solutions of the four-dimensional theory correspond to instantons of the as- 
sociated three-dimensional non-linear sigma model over G5/fj*. The spherically symmetric 
black holes (including the asymptotically Taub-NUT ones) are entirely characterised by 
their (5-Noether charge 

Q = -L / +VPV- 1 (2.4) 

and the asymptotic value of the scalars fields V„ G 64 C at spatial infinity [6] . Here 
it will be more convenient to characterise the solutions in term of a modified conserved 
charge ^ obtained by rotating Q back into the coset 

tf = V - l QV eget)* (2.5) 

which we will call the 'Noether charge' for simplicity (this designation being unambiguous 
since we will never refer to Q itself). 

Any spherically symmetric black hole then admits the following representative 

[19] 

v =Ws in (^H (2 - 6) 
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for some matrix where r is the Weyl radius, associated to the three-dimensional 
Riemannian metric 1 

g^dafdaf = dr 2 + (r 2 - c 2 ) {d6 2 + sin 2 Bdip 2 ) (2.7) 

The BPS parameter c is defined in terms of the Noether charge as 

c 2 = \ Tr ^ 2 , (2.8) 

Here k = Tr h 2 is a positive integer associated to the group © and the following five- 
graded decomposition of its Lie algebra g with respect to its subalgebra gt x = Rh 

1^ © © © f ) © © 1( 2 > (2.9) 

This decomposition shows explicitly how the duality group (& 4 of the four- dimensional 
theory and its representation t 4 are embedded in (5. 

Given a solution of the three-dimensional equations of motion one can reconstruct the 
solution in four dimensions as follows. In order to read off the Kaluza-Klein ansatz for 
the four-dimensional metric and the abelian vector fields from the above 'data', one must 
first rotate the coset representative V into a parabolic (triangular) gauge V G ^3 where 
*P is the maximal parabolic subgroup of & defined from the five-graded decomposition 
(2.9) with Lie algebra p 

P = (fltiffiflf) (2.10) 

Let us designate by $ the involution that defines ()* C g, such that x — x* G f)*. We write 
h for the generator, and e and f = e* for the grade 2 and —2 generators, respectively. 
The representation [ 4 always admits a (55 4 invariant symplectic form (which defines the 
Dirac quantisation condition) such that two generators x and y of [ 4 commute to 

[x,y] = (x,y) u e (2.11) 

Defining $ as the electromagnetic scalars valued in 1 4 and v G 4 /i} 4 the four-dimensional 
scalars fields, one writes down the coset representative V in the parabolic gauge as 

V = exp( J Be + $) expQln/fhJ v (2.12) 

lr The Weyl radius r is related to the standard radial coordinate f via f = r + m such that e.g. the 
horizon of the Schwarzschild solution corresponds to r = to (sec e.g. [20]). Unlike the usual Schwarzschild 
coordinates the Weyl coordinates are duality invariant, as follows from the duality invariance of the BPS 
parameter c and the split of the four-metric in (2.15) below. 
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A straightforward computation gives 
P = l -H- l dH h + l -H- 1 (dB - l - ($, d$) J (e + f ) 

+ itf-il^tr 1 )^* + **) + (tT 1 *;) |04Qh4 (2.13) 

where L 4 is the l 4 representation homomorphism. Substituting these expressions into the 
equation of motion (2.2) it follows that the axion field B defines the Kaluza-Klein vector 
B through 

H- 2 * [dB- -($,d$)J = dB (2.14) 
such that the four-dimensional metric reads 

ds 2 = -H(dt + B^dx*) 2 + H- 1 g^dx^dx" (2.15) 

In order to reconstruct the four- dimensional vector fields 

V2ttGA = U(dt + B^dx* 1 ) + (2.16) 

from (2.13) one must choose a Lagrangian subspace of [4 with respect to the symplectic 
form of <3 4 , such that $ splits into U © A. After the redefinition B' = B + |(?7, A)^, the 
fields A only enter P linearly through their differential dA, and can then be dualised to 
vectors A according to their equations of motion. 

Let us consider a general non-rotating asymptotically flat solution which does not 
carry any naked singularity, by which we mean that any singularity is covered by an 
horizon which is identified with its corresponding Killing horizon (no ergosphere). Such 
solutions are spherically symmetric black holes [6] . From the point of view of the three- 
dimensional Riemannian space V, the horizon J4? corresponds to the single point r = c on 
V (cf. (2.7)) with an instanton-like singularity. Because the space-time volume element 
H^^fg is regular on the horizon of a non-extremal black hole, while H goes to zero 
there, the three-dimensional volume element ^fg vanishes on the horizon. It follows that 
all components of V must tend to the same value as r — > c and are thus constant on M' . 
Using the fact that the dependency on B can be neglected on 3% ', the expression of the 
horizon Aj? in the coordinates (2.7) 

A^ = I (r 2 - c^H' 1 siuddd Ad<p (2.17) 
J.je 

shows that for Aj? > 0, H behaves like 

H = i^(r 2 - c 2 ) + 0((r 2 - c 2 )l) (2.18) 
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near M' . Similarly, the surface gravity K can be computed from the Killing vector 
£ M d M = d t as 



K 2 = -^im V m £ n Vm£n = j 



dWdfJI - H d^B v {d^By - d v B») (2.19) 



Exploiting the fact that B is regular on the horizon we thus obtain 

k = -d~H 
2 



V= — V 
An 



(2.20) 

\je Ajjf 

Similarly, using the behaviour of the fields near the horizon, one computes that the charge 
associated to the horizon is given by 

— Vo' 1 [ *VPV~ l V = lim ( (r 2 - c^V^VPrV^V^ (2.21) 

It follows that 

Tr <^ 2 = lim ({r 2 - c 2 ) 2 Tr P 2 ^j (2.22) 

Using the behaviour of H near the horizon, the equation of motion (2.14) and the as- 
sumption that all the fields B, $ and v are regular on the horizon, we arrive at 

(r 2 - c 2 ) P r = c h + O (Vr 2 - c 2 ) (2.23) 

recovering the formula (2.8) as expected. Therefore, we obtain that the product of the 
horizon area and the surface gravity of a non-rotating black hole is determined by the 
Killing norm of the 0-Noether charge as 

A^k = 4ttc (2.24) 

As a result, non- rotating extremal solutions, for which K = 0, carry a Noether charge 
satisfying Tr ( € 2 = (since K is identified with the temperature of the black hole in the 
thermodynamic interpretation, it follows that all such solutions are characterized by zero 
temperature). Therefore, from (2.7) we see that spherically symmetric extremal black 
holes solutions can be viewed as instantons over flat three-dimensional space, with the 
Euclidean three-dimensional metric g^ u = 5^ u . For these solutions the coset representative 
V takes the form 

V = V,exp(~ < iA (2.25) 

with a Noether charge satisfying Tr c & 2 = 0. Inspection of the equations of motion 
shows that - can be replaced by any solution of the three-dimensional Laplace equation 
AH. = 0. Such solutions were first discovered by Papapetrou and Majumdar in Maxwell- 
Einstein theory [1, 2]. The main purpose of this paper is to study these solutions in a 
systematic manner by exploiting the maximal duality symmetry & acting on them. 
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3 Solutions of Papapetrou— Majumdar type 



To obtain the solutions of Papapetrou-Majumdar type, we replace the formula (2.25) by 
the general Ansatz for the coset representative V 

V(x)=V exp (-^2n n {x)^A (3.1) 

with Lie algebra elements ^ £ g0 f)*, and where the functions 7i n (x) are to be determined 
from the equations of motion. It is straightforward to see that P depends linearly on the 
functions 7i n if and only if [9] 

[V m ,[V n ,%]]=0 (3.2) 
Under this assumption, one obtains that 

P = ^dU nc € n 

and the equations of motion reduce to 

= d^H'd v W Tr d * dH" = (3.4) 

mn 

For the equations of motion to be linear in Ti" , one must require in addition that 

Tr = (3.5) 

Any set of matrices satisfying (3.2) and (3.5) yields a solution of the theory with a 
flat three-dimensional metric, if the 7i n are arbitrary harmonic functions. Therefore the 
problem of solving the three-dimensional field equations can be reduced to solving the 
algebraic equations (3.2) and (3.5). The main new insight from the present analysis (and 
from [8]) is that this problem, in turn, can be reduced to the construction and classification 
of the nilpotent orbits of q via their corresponding normal triplets (see below) and the 
associated graded decompositions of q, as we shall now explain. 

Any element of a Lie algebra g can be written as the sum of a diagonalisable element 
and a nilpotent element. There is thus no loss of generality in assuming that each 
c € n in the general ansatz (3.1) is either nilpotent or diagonalisable. Two elements ^ 
and ^ satisfying equation (3.2) either commute, or generate a Heisenberg subalgebra 
of g. From the representation theory of Heisenberg algebras it follows that any linear 
combination of c € 1 and ^ is nilpotent in the latter case. Therefore the matrices c € n are 



Q = - ^2 'WdW 1 [% , c € m 



(3.3) 
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either diagonalisable elements commuting with all the others, or nilpotent elements such 
that any linear combination of them is likewise nilpotent. 

The set of nilpotent elements of g is known to decompose into finitely many ©-orbits, 
the so-called nilpotent orbits of g [15]. For a given nilpotent orbit, one can associate to 
a given representative E G g of the orbit, an triplet (H, E, F) such that [15] 

[H, E] = 2E [H, F] = -2F [E, F] = H (3.6) 

Decomposing the Lie algebra with respect with the eigenvalues of H, one gets a (2n+ 1)- 
graded decomposition of g (where g[ x = KH) 

g <-«) © ... © fl <- 2 > © gt" 1 ' © (g^ © g<°>) © g« © g< 2 > © ■ ■ ■ © S '"» (3.7) 

such that g (0) is a reductive Lie algebra (i.e. the direct sum of a semi-simple and an 
abelian algebra) which contains a Cartan subalgebra of g. From equation (3.6), it follows 
that the nilpotent element E lies in g (2) , and the nilpotency degree of E in the adjoint 
representation 2 is the lowest integer strictly greater than n. The nilpotency degree then 
follows directly from the grading and is not altered by the detailed commutation rules 
of the Lie algebra in this graded decomposition. For a general nilpotent orbit, n can be 
pretty large, but if we consider a nilpotent orbits for which n < 5, it follows from the 
graded decomposition that any element ^ G ©™ =2 <p) satisfies both equations (3.2) and 
(3.5). There are four distinguished cases, namely 

ad E 6 = n = 5 

5 



ad E =0 n = 4 

1 E 4 



ad E 4 = =>• n = 3 



ad E 3 = n = 2 (3.8) 

All of the above relations also hold for the complexification of the duality group © and 
its associated complex Lie algebra Qc- What is important is that the complex nilpotent 
orbits in gc are uniquely determined by the nilpotency degree (in several representations), 
and equivalently by the graded decomposition (3.7). When descending from the complex 
to the real Lie algebra, the complex nilpotent orbit may decompose into several real 
orbits, which are the ones relevant for black hole solutions. 

The isotropy subalgebra of a representative E of a given orbit is a subalgebra of 
0™ =o g (p) . Any solution to equation (3.2) can thus be associated to a given 2n + 1 graded 



2 Which is defined to be the smallest integer k such that (adE) fc (x) = for all x £ g. 
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decomposition (3.7) with n < 5, such that the nilpotents elements lie in ©™ =2 g <p) and 
the diagonalisable elements commute with all the others and lye in g (0) . Equation (3.5) 
is trivially satisfied for the nilpotent elements, but strongly constrains the diagonalisable 
elements of g (0) . 

More specifically, in black hole physics, we are interested in elements which lie in 
g h* in order for the coset representative V to be in the symmetric gauge. In fact, one 
can show that there is no loss of generality by doing so, because any more general coset 
element can be rotated back to this form by a right fy* gauge transformation. When the 
(J5-orbit of the nilpotent element E has a non-trivial intersection with g h*, the triplet 
can be chosen in such a way that both E and F lie in g ()*, and such that H lies in f)*. 
It follows then that h* decomposes in a similar way as 

(f |j<-> © ... © h ( - 2) © fj^ © (g[ x © h (0) ) © h (1) © h (2) © ••• © f)<"> (3.9) 

such that f) (0> is a reductive Lie algebra which contains a Cartan subalgebra of h*, and 
E G g (2) f) (2) . The elements of g (p) h (p) (for p > 2) generate a nilpotent Lie algebra 
n (p) , which decomposes as (for n < 5) 

n (p) = (g (p) © h (p) ) © h (2p) (3.10) 

Any set of charge matrices ^ lying in ®" =2 (g (p) 0f) (p) ) therefore satisfies the commutation 
relations (4.18) and defines a multi-black hole solution via the Ansatz (3.1). In the case 
of the graded decomposition (2.9), which will be the one of interest as we are going to 
see, the isotropy subgroups of grade zero 3l 0> C & and ^ 0> C $)* of a generic element 
E G [4 f) (2) coincide with the isotropy subgroups of <&4 and ^4 of the corresponding 
electromagnetic charges, respectively. These subgroups also define the so-called moduli 
spaces of black hole attractors as ZTI^-T c ^4/^4 [14]. The Cartan norm is thus 
positive definite on j^ 0) and there is no diagonalisable element to consider in the 
solutions of (3.2, 3.5). 

Let us next consider a more general Ansatz than (3.1) including Lie algebra elements 
■s^s G ©p =2 f) (n) - From the grading and the Campbell-Hausdorff formula it follows that 

V = V exp (- ^ n"^ n - K s srf\ 

= V exp(-Y^n*K + lj2 VH *^^ exp^-^/CW^ (3.11) 
such that this Ansatz is in fact equivalent to (3.1) up to a right $)* gauge transformation. 3 



3 Note that kx] e e; =4 (0 (p) © f) <P) )- 
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The solutions obtained by the above construction may in principle exhibit naked 
singularities. Quite generally, the latter can be of two types. For solutions not obeying 
any 'no force' property, there are usually rod-like singularities which appear in static 
multi-black hole solutions in order to balance the gravitational pull between any two 
black holes, and which manifest themselves as string-like singularities of the Riemannian 
metric g^. For example, in the axisymmetric case, 4 they appear in Weyl coordinates 

g^dxPdx" = e 2a (dz 2 + dp 2 ) + p 2 d^ 2 (3.12) 

if the function a(p, z) does not vanish in the limit p — > on the axis in between two 
interacting black holes. This type of singularity cannot occur for the solutions built on 
the conditions (3.2) and (3.5) because the three-dimensional metric is then flat and regular 
everywhere by (2.3). Consequently, the singularities of the four- dimensional theory solely 
originate from the singularities of the scalar fields which are located at the poles of the 
harmonic functions H,". Naked singularities of such solutions are thus entirely due to 
the individual black holes themselves, and can be avoided by choosing each charge c € n 
so that the individual black holes have their singularities covered by horizons. This way 
one makes sure that each pole of the harmonic functions corresponds to an horizon in 
four dimensions. 

It is commonly assumed that all the spherically symmetric extremal black holes with- 
out naked singularities correspond to particular limits of non-extremal spherically sym- 
metric black holes which do not carry any naked singularity themselves. Within the 
class of model discussed in this paper, the asymptotically flat regular non-extremal so- 
lutions are all in the fj*-orbit of uncharged Kerr-solutions with scalar fields having any 
constant value [6]. It follows that the corresponding Noether charge ^ satisfies a cubic 
characteristic equation [8] 

^ 3 = c 2 ^ (3.13) 

in the fundamental representation (which is the spinor representation of its double-cover 
when & is an orthogonal group), save for two particular cases for which (25 is a non- 
compact real form of E$, corresponding to maximal M = 8 supergravity and the 'magic' 
M = 2 supergravity associated to the octonions [21]. In these two cases the characteristic 
equation is quintic [8], and must be satisfied in the 3875 representation that appears in 
the symmetric product of two copies of the adjoint representation, 

^ 5 = 5cV 3 -4c 4 ^ (3.14) 

In which case exact solutions with several black holes that do not share any no force property are 
known explicitly [20]. 
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It follows that regular extremal solutions (for which c = 0) carry a nilpotent Noether 
charge which vanishes at the third power in the fundamental representation, i.e. 
c )o 3 = 0, in the generic case, or at the fifth power in the 3875 representation of E$ when 
is a real form of E 8 [8]. So as a first consequence, one can restrict attention to nilpotent 
elements in considering regular multi-black holes, because any diagonalisable elements 
of vanishing Killing norm would lead to naked singularities. 

The nilpotency conditions are preserved by the action of the complexified group 0c, 
and they determine one single complex nilpotent orbit with representative E 5 . As already 
mentioned the complex orbit may decompose into several real nilpotent orbits when one 
descends to the real group C 0c- Acting with the complex Lie group 0c on the 
representative element E 5 we obtain a dense open subset within the complex variety 
defined by the nilpotency condition [15]. Therefore, any solution ^ of the characteristic 
equation (3.13) or (3.14) lies in the closure of this complex orbit inside the variety of 
nilpotent elements of 0c- Since, in addition we require this orbit to be real and to have 
a non-trivial intersection with Q Q f)*, we have 

^e6 c -E 5 n ef)* (3.15) 

It turns out that the sl 2 (C) triplet associated to E 5 is such that H; = [E J7 F 5 ] satisfies the 
same characteristic equation as 2h, and the graded decomposition of 0c associated to such 
nilpotent element is therefore the (complexified) five-graded decomposition associated to 
the dimensional reduction (2.9). The five-graded decomposition of 0c with respect to H 5 
corresponds to a unique five-graded decomposition of 0, 5 

g T- 4 > © l[~ 2) © (0[ x © f ) © [f © l« (3.16) 

In general, the elements E and H of a normal triplet do not satisfy the same characteristic 
equation, and this property is very particular to h and its five-graded decomposition. 
In fact, this five-graded decomposition characterises both the minimal semi-simple 0c- 
orbits of 0c, i-e. 0c ■ h — 0c/(C x x 04c), and the minimal nilpotent c -orbit of 0c, 
i.e. c -e = 0c/(04cx(ti 1 c ©C^)). 

In order to correspond to a regular solution, the Noether charge must lie inside the 
closure of the ij*-orbits inside 0c • E 5 fl0 ()*, whose isotropy subgroups are contracted 
forms of #4 [8], since the latter is the isotropy subgroup of the Kerr solutions. S)^ being 
compact, the semi-simple component of the isotropy subgroup of Sj* of such charges is 

5 Although it can be degenerate, as for example for e 6 (_i4) and f7(_25)- Nevertheless the orbit we are 
interested in is unique in these cases. 
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also compact. The most general solutions of Papapetrou-Majumdar type correspond to 
linear combination of such nilpotent charges which vanish at the sixth power in the adjoint 
representation. When considering a nilpotent element within the graded decomposition 
associated to a higher order orbit, its isotropy subgroup generally involves components of 
negative grading such that the semi-simple component of its isotropy subgroup appears 
itself through a graded decomposition embedded inside the graded decomposition of h* 
associated to the higher order orbit. It follows that the latter isotropy subgroup cannot 
be a contracted form of ft 4. The relevant linear combinations of charges therefore all lye 
in the closure of the Sj* orbits of 0c • E 5 D g h*. We will prove explicitly in [13] that 
this is indeed the case within M = 4 supergravity coupled to n vector multiplets. The 
only nilpotent orbits we have to consider in order to obtain the most general solutions of 
Papapetrou-Majumdar type are thus the real nilpotent orbits of g of dimension 2 dim^]-!- 
2 associated to the five-graded decomposition (3.16), or more precisely, the dim [[4] + 1 
dimensional f)*-orbits of nilpotent elements inside g h* associated to this five-graded 
decomposition, such that the subalgebra lies in ()*. In general, there are different 
inequivalent embeddings of f)* C g consistent with the graded decomposition (3.16), 

f>* f,C- 4 ) © f,C" 2 ) © glj © f,<°> © |j< a > © J,« (3.17) 

which correspond to the different i}*-orbits inside 0c ■ E 5 H g f)*. There are inequiva- 
lent Papapetrou-Majumdar solutions associated to each of those graded decompositions, 
whose number of independent harmonic functions are dim^ h (2) ] + 1 — dim[h (4) ]. 



4 Supersymmetric solutions 

When the theory coincides with the bosonic sector of a supergravity theory, it is well 
known that BPS black holes preserving identical supersymmetry charges define multi- 
particle solutions of Papapetrou-Majumdar type. We are now going to see how this 
can be understood in the general framework developed in the foregoing section. To this 
aim we recall a main result of [8], namely the fact that the Noether charge admits a 
dual description, either as a matrix c io valued in g f)* as in (2.5), or alternatively, 
and equivalently, as a state |^) in a fermionic Fock space, transforming in the same 
representation as the coset matrix. In the latter description the BPS parameter (2.8) is 
simply given by the indefinite 'norm' of this state, to wit 

c 2 = (^|%?> (4.1) 
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For TV-extended super gravity, the non-compact group Sj* is the product of the chi- 
ral component 6 of the R-symmetry group Spin*(2J\f) and a group SJq depending on 
the matter content of the theory. The supersymmetry parameters are SU (2) spinors 
valued in the pseudo-real vector representation of SO*(2Af) which satisfy an SO*(2Af)- 
Majorana condition associated to the pseudo-anti-involutions of SU(2) and SO*(2Af). 
Using U(N) C SO*(2J\f) covariant notations, the spinor parameters are complex SU(2) 
spinors e l a valued in the fundamental representation of U(Af), with complex conjugates 
ef. In order for the solution to be supersymmetric, the corresponding Noether charge 
state must satisfy the 'Dirac equation' [8] 

(e^ + e^efa*)!^ = (4.2) 

where a 1 and a, (for ■ ■ • — 1, ...,7V) are the fermionic oscillators from which the spinor 
representations of SO*(2J\f) are built (see appendix B of [8] for details). The condition 
(4.2) is derived by considering the supersymmetry variation of the dilatino fields in the 
asymptotic region. Observe that it contains more detailed information than the gravitino 
variation, which simply yields the necessary condition c = for Killing spinors to exist 
(the latter statement also follows from (4.2) as an integrability condition). In searching 
for extremal solutions, we can thus replace the problem of solving the cubic or quintic 
nilpotency conditions for the charge matrix ^ by the simpler linear condition (4.2) for 
the associated charge state \€) [8]. 

Solutions of (4.2) are characterised by the number of supersymmetries and by their 
isotropy groups. Generally, the latter are always parabolic subgroups of Sj* because the 
charge state |^) is a null vector (i.e. a zero norm state) for c = by (4.1). In [8] we 
showed that solving the characteristic equations (3.13) or (3.14) for regular solutions 
gives rise to a stratified moduli space 

M = MoUM^-'-UMn (4.3) 

where the main stratum Aio corresponds to non-extremal solutions (that is, c 2 > 0), while 
the remaining strata correspond to solutions with c = and contain the supersymmetric 
or BPS solutions as a subspace. When the solution is left invariant by 4n supersymmetry 
generators, one speaks of an BPS solution'. In this case a basis of supersymmetry 
generators can be chosen with spinor parameters with A running from 1 to 2n, which 

6 By chiral component of Spin* (2J\f), we mean its subgroup Spin*(2J\f)/ kcr(S' + ) that acts faith- 
fully in the chiral Weyl spinor representation S+, i.e. SL(2,WL) C Spin*(4), ,50(2,6) C Spin*(8), 
Spin* (16) /Z 2 C Spin* (16) and Spin*(2Af) itself for odd N. 
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satisfy the reality condition (alias symplectic Majorana condition) 

e A + £a pn AB ei = (4.4) 

for a given symplectic form Q AB of C 2n , satisfying Q A c^l BC — $a- Using (4.4) the 'Dirac 
equation' (4.2) can be further simplified to 

(a A -n AB a B )\^) = (4.5) 

Let us designate by ?P« the parabolic subgroup of SO*(2Af) that preserves the set 
of solutions of (4.4), and thus the set of supersymmetry charges leaving invariant a 
solution of Noether charge \ c io) obeying (4.5). To derive we note that the Lie 

algebra associated with *P « is generated by the eigen vectors of positive eigen value with 
respect to a g^ generator D« satisfying D^.e A = e A , which reads 

Dn. = - (Q AB a A a B - Q AB a A a B ) G so*(2Af) (4.6) 

in the harmonic oscillator basis (see [8] for details). Its action on the spinor parameters 
satisfying (4.4) can be read off from 

[ Dn, , e A a A + e af3 e p A a A ] = e A a A + e aP e p A a A (4.7) 

D» defines a five-graded decomposition of S0*(2JV), 

S0*{2M) = n(2n - 1) (_2) © (2n <g> 2n ± ) ( R 1) 

© (g[ x © su*(2n) © so*(2n ± )) (0) © (2n © 2n ± )« © n(2n - 1) (2) (4.8) 

with n± = M — 2n. Indeed, using the identity, 

[D j L 7 a A ± n AB a B } = ± (a A ± Q AB a B ) (4.9) 

one concludes that the grade zero subalgebra is generated by 

^A A B [a A + {l Ac a c , a B + Q BD a D ] G su*(2n) (0) 

U A B [a A ,a B ] + h AB a A a B - h AB a A a B G *o*(2n x )<» (4.10) 

where A is a U(n±) index running from 2n + l to TV. The nilpotent generators in so* (27V") 
of grade ±1 and grade ±2 are, respectively, 

A AB (a A ±n Ac a c )a B - A AB (a A Ttt Ac a c )a B G (2n © 2n ± )^ 1) 
l -A AB (a A ± n Ac a c ) (a B ± f^ D a D ) -V s (a A T n Ac a c ) (a B T n BD a D ) 

G n(2n-l) (±2) (4.11) 
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The covariance of the spinor parameter satisfying (4.4) under 5u*(2n) follows from the 
commutation relations 

{a A + Q AC a c , a B + Q BD a D } =0 {a A + Q Ac a c , a B + Vt BD a D ) = 28 A (4.12) 

The parabolic subgroup C SO*(2Af) preserving the set of spinor parameters satis- 
fying (4.4) is thus defined by the Lie algebra 

p% = (sk © su*(2n) © so*(2n ± )) (0) © (2n ® 2n J _)£ ) © n(2n - 1) (2) C so*(2Af) (4.13) 

This group preserves the set of spinor parameters (4.4) in the sense that they are only 
rotated by the subgroup SU*(2n) and rescaled by GL(1,R), while the remaining gen- 
erators of act trivially on them. In conclusion, acting on an jj BPS spherically 
symmetric solution with an element of ^3 « x S)q, one obtains another BPS spherically 
symmetric solution which preserves the same supersymmetry charges. 

For asymptotically flat solutions, i.e. asymptotically Taub-NUT space-times, the <3 
Noether charge ^ is valued in 00f)*. In supergravity, this is an irreducible representation 
of Spin*(2J\f) xS)q, that is, the product of the Weyl spinor representation S + of Spin*(2Af) 
and an irreducible representation TZ of $Jq, possibly with an extra reality condition. 

For simplicity let us first restrict ourselves to -k BPS solutions. The theories admit- 
ting BPS solutions of higher BPS degree are quite restricted anyway and further relevant 
examples will be discussed in [13]. In this case the generator defines the five graded 
decomposition (which is an example of the decomposition (3.17)) 

so*(2A/") = 1 ( ~ 2) © (2 © 2n ± ) i k 1) © (flli © su(2) © so*(2n ± )) (0) © (2 © 2ni)£ © 1 (2) (4.14) 

and the Weyl spinor representation of Spin*(2Af) decomposes accordingly as 

S + S^ 1] © (2 © S_) m © (4.15) 

where S± are the Weyl spinor representations of SO*(2n±). Using the fact that the 
-^=-BPS Killing spinors lie in the grade-one component of the vector representation of 
SO*(2Af), 

2N = 2 ( - l) © 2n ± (0) © 2 (1) (4.16) 

one sees that the solutions of the Dirac equations are simply obtained by taking the 
charge matrix ^ to lie in the grade-one component of g f)*, 

(,S + ©ft) (1) CS+®K (4.17) 
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It follows from the graded decomposition of g with respect to Dj_ that for any charge 
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matrices % G (S+ ® ft) (1) , 

[tf n ,K] e 1 (2) C so*(2A0 =*► [<jf„ = ( 4 - 18 ) 

In this way one obtains that such solution is of the type discussed in the preceding section, 
with H = 2D i . 

77 

For n > 1, the general solution of (4.5) is obtained by use of the identity 

(a A - n AB a B ) exp f -Sl CD a c a D ) |0) = (4.19) 



2 

Consequently, if the charge ^ is represented by the state 



\<#) =[(w + Z i3 a l a? + • • •) |0) , [z A + E^aV + • • •) |0) ^ (4.20) 



the general -j^-BPS solution reads 



|#) = ^ e * nABaAaB (w + Z AB a A a^ |0) , e ^ ABaAaB z A |0) j (4.21) 

where A runs from 2n + 1 to A/", and >1 labels the matter multiplets if M < 4. An 
additional reality condition (Majorana-Weyl condition) may have to be imposed on the 
state if 00f)* admits one. It ensures that the scalar charges are completely determined in 
function of the others for any solution of the Dirac equation. 7 In this case, the generator 
H is identified as 

H = ~ D ^r ■ (4-22) 
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and E with the charge matrix ^ in such a way that the commutator [H, E] corresponds to 
the action of the operator (4.22) on the state |^). 8 Because the component of grade two 
(with respect to D^) of so*(2A/*) leaves the jJ=-BPS Killing spinors invariant, the Killing 
equation trivially reduces to de A = 0, and it follows from (4.5) that the corresponding 
supersymmetry variation of the dilatino fields vanishes. 

Having constructed the appropriate charge matrices ^ it is now straightforward to 
construct multi-black hole solutions of the four- dimensional theory following the proce- 
dure described in section 2. If the different matrices c io m satisfy (4.2) with respect to 



7 For instance for M = 2, 3 the scalar charges are manifestly determined in (4.21). For AT = 4 the 
complex self-duality condition on the vector multiplets ensures that this is still the case, and similarly 
for J\f = 8 through the complex-selfduality of the multiplet. 

8 That \ c (o) is always null for BPS states obeying (4.5) then simply follows from Tr E 2 = 0. 
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the same supersymmetry parameters, the resulting multi-black hole solution will also 
be supersymmetric. If the 'constituent' black holes are supersymmetric with respect 
to different supersymmetries, on the other hand, the multi-black hole configuration will 
no longer be supersymmetric, and in general then exhibit naked singularities unless the 
'constituent' BPS black holes have a vanishing horizon area. Examples of such solu- 
tions which are regular are given in [13]. Although supersymmetry is not required for 
the existence of Papapetrou-Majumdar type solutions, it greatly simplifies their explicit 
derivation. This is because the algebraic equations determining representatives of nilpo- 
tent orbits are not trivial to solve in general, whereas the linear Dirac equation (4.5) is. 
Nevertheless we will see that the supersymmetric graded decompositions can be used to 
derive the explicit form of the nilpotent charges in the non-supersymmetric well. 

5 Non-BPS solutions in five- dimensional theories 

The second class of solutions appear in the theories that can be obtained by toroidal com- 
pactification of (super)gravity theories in five dimensions to four dimensions. Inspection 
of the list of theories displayed in [6] shows that this is the case whenever $)* is semi- 
simple. The scalar fields of the five-dimensional theory then parametrise a symmetric 
space ®5/f)5, and the isometry group ©5 acts on the abelian vector fields in a faithful 
representation l 5 . One then has the following three-graded decomposition of 04 and the 
associated four-graded decomposition of [4, 

S4 [<r 2 > © © gf ) © [f [ 4 It"') © [<-) © © 1<3) (5.1) 

The maximal compact subalgeba fj 4 = f) 5 © [ 5 admits a non-compact real form f) 4 = 
f)5 ® U5 which appears in the pseudo-Riemmanian coset space <8a/$)1 of the Euclidean 
four- dimensional theory obtained by dimensional reduction along the time-coordinate. 
Using the graded decomposition one obtains the following result: because l 4 defines a 
complex irreducible representation of f) 4 C 34, it decomposes into two conjugate real 
irreducible representations, l 4 and T 4 , of f) 4 C g4. This property can be understood as 
follows. Let us define the four- dimensional field strength T in the 1 © [5 representation, as 
well as Q = 5C/ bT (which in general depends on the scalar fields) such that the equations 
of motion of the abelian vector fields are d-k Q = 0. The fj 4 -invariant form of the action 
on a four-dimensional pseudo-Riemannian space-time involves the complex combination 
*G + iT which transforms in the complex representation [4 of ^4. This representation 
is complex because * 2 = — 1 for a pseudo-Riemanian metric, whereas the fj 4 -invariant 
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form of the action on a four-dimensional Riemannian space involves the two independent 
real combinations *Q ± T which transform in the real representation [4 and its dual [4, 
respectively. As an example consider M = 8 supergravity, where -kQ + %T belongs to 
the complex representation 28 of SU(8), whereas in the Euclidean signature theory the 
corresponding combinations -kQ ± T belong to conjugate real representations 28 and 28 
of SU*{8). 

It follows that when the four- dimensional theory descends from a five- dimensional one, 
the Lie algebra f)* — 5I2 © f)s © 3 <g) [5 admits the following three-graded decomposition 

f,* 9a*c- 2 >©( [ lffi I)*)C°>©"[f) (5.2) 

which gives an example of decomposition (3.17), with t) (4) = in this case. Indeed one 
obtains the five-graded decomposition of g © f)*, 

e r = i ( ~ 4) © ^T 2) © (04 © f)D (0) © iT © i <4) (5-3) 

such that I4 © 1 (4) defines an abelian nilpotent subalgebra of g. One can thus define 
solutions of Papapetrou-Majumdar depending on dim[[ 5 ] + 2 independent harmonic func- 
tions. As we will see, the most general asymptotically Minkowski solutions of this type 
are defined in function of charges lying in [4 and depend on dim [[5] + 1 independent 
harmonic functions, and the extra-harmonic function corresponds to considering these 
extremal black holes in the background of pure NUT extremal black holes with vanishing 
horizon area. 

6 Maximal supergravity 

We now wish to illustrate the general discussion of the last sections through the exam- 
ple of the maximally supersymmetric supergravity. This theory being also the toroidal 
dimensional reduction of the maximally supersymmetric supergravity theory in five di- 
mensions, it carries the two kinds of Papapetrou-Majumdar solutions discussed in the 
two preceding sections. They define the most general solutions of Papapetrou-Majumdar 
of the theory, which correspond to the only two £Vs) nilpotent orbits associated to the 
graded decomposition (3.16), respectively [8, 23]. 

Let us first recall the form of the E&rg\ Noether charge in the harmonic oscillator basis 
[8] , with the mass and NUT charges w = m + in, the electromagnetic central charges Zij 
and the scalar charges The latter are determined as functions of w and through 
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the characteristic equation (3.14). The state \€) G Cs(8) 0so*(16) is explicitly given by 

= (w + Zjja'a 1 + Z ijk ia i a j a k a l + ^e ijk i mnpq Z pq a { ■ ■ ■ a n + ^e ijk i mnpq w a? ■ ■ ■ a q ^j |0) 

= (l+*)(w+ %aV + -EyjuaVW) |0) (6.1) 

where ~k is the anti-involution defining the 12 8- dimensional chiral Majorana-Weyl repre- 
sentations of Spin* (16) [8]. In the following subsections we will exhibit the form of \ c &) 
explicitly for various kinds of extremal black holes. 

6.1 \ BPS solutions 

For a generic | BPS solution preserving the supersymmetry generators associated to the 
Killing spinor parameters 

4 + = 4 = ( 6 - 2 ) 
with A running from 3 to 8, the Noether charge must satisfy 

(a 1 + a 2 ) \<g) = (a 2 - a x ) \<g) = (6.3) 

The general solution can be written as 

= (1 + *) (l + aV) (w + Z AB a A a B ) |0) (6.4) 

This spinor evidently obeys the Majorana-Weyl condition ~k\io) = Y^), and to see that 
it still satisfies (6.3) in accordance with (4.19) one uses the identity 

(1 + *) (l + aV) (w + Z AB a A a S ) |0) = (l + aV) (1 + *) (w + Z AB a A a B ^j |0) (6.5) 

where -k is the anti-involution for so*(12). The grading generator Di G so*(16) is 

Di = a}a 2 — dido (6.6) 

8 

and obeys 

D, \V) = |tf> (6.7) 
for any value of Z AB . It gives rise to the following five-graded decomposition of so*(16) 

so* (16) = l ( - 2) © (2 © 12)^' © (q{ x © su(2) ©so*(12)) (0) © (2 © 12)™ © 1 (2) (6.8) 

The associated charge ^ G ts(8) ©so* (16) is left invariant by su(2) and transforms as a 
Majorana-Weyl spinor 32 + of so*(12), with the decomposition 

e 8(8 ) © so* {Id) = 32^ 1) © (2 © 32_)£ 3 © 32™ (6.9) 
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To see this explicitly, we invoke (6.5) once more. The action of the grade zero component 
R*j_ x SU(2) x Spin* (12) C Spin* (16), preserves the given supersymmetry charges by 
construction. In a similar way as for the action of Spin* (16) on the Kerr solutions, the 
moduli space of regular black holes preserving these four supersymmetry charges is the 
closure of the R+ x Spin* (12) orbit of a given regular solution that can be chosen to 
correspond to the charge matrix = (l + a 1 a 2 )(l + *)|0). This orbit is 9 

, Rix Spin*(12) , A Spin*(16) 

suW ~ cM ^ HsmhuW) (6 ' 10) 

This embedding is a refinement of the stratification (4.3), with the following notation: 
by M.™ for m < n we generally designate the space of -j^-BPS solutions which are 
such that m out of n supersymmetry charges are shared between all of its elements 
(which implies M.™ C JA™ for m > m'). In other words, while A4± contains all |-BPS 
solutions, irrespective of which supersymmetry charges are left invariant, A4\ consists of 
all solutions that are | supersymmetric with respect to the same supersymmetry charge. 
Consequently, superposing solutions with charges c € m from the orbit M\ according to 
the basic ansatz (3.1) will result in a multi-black hole solution that still preserves |- 
super symmetry. By contrast, taking different c € m from the larger orbit Aii in general will 
lead to a solution that is no longer supersymmetric because the | supersymmetries of its 
black hole 'constituents' do not match. 

Positivity of the energy entails that \w\ is strictly greater than any of the eigenvalues 
of Zab, an d that the quartic £7(7) invariant associated to it is strictly positive. Note that 
to avoid naked singularities, one must require moreover w G C \ R_, and w G for 
asymptotically Minkowski solutions. The boundary of the orbit (6.10) decomposes into 



9 See [8] for our notations (which are not standard): by '7c' we designate the parabolic extension of 
SU(2) x SU(6) by the grade-one translation generators and the grade-two central charge 

Ic(SU(2) x SU(6)) = (SU(2) x SU(6)) x ((2 <g> 6) (1) © 1 (2) ) 

Note also that the SU(2) factors in the quotient defining M\ in (6.10) 'cancel out'. 
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orbits corresponding to black holes with vanishing horizon area 10 

Spin* (12) /- KA ~ Spin* (16) 

jvl l° — Sp(3)xl4 2 ^jvi la — R * Kj c (ST/(2)xSp(3)tKl4 2 ) 

A// 1 ~ Spm*(12) Spin* (16) ((\ ~\'\\ 

JV1 2 — (5C/(2)x5pm(l,6))K((2<g)8) R ©l) ^ J 12 — (5J7* (4) xSpm(l,6)) K((4®8) R ©6) V u - xj V 

A^fl ~ Spm*(12) f' — A/f — Spin* (16) 

yVl 4 — SC/*(6)xl5 ^ ;vl4 — SC/*(8)x28 

with the notation from [8] (see especially eq. (5.17) there) for the embedding of the various 
'BPS strata' into the stratified space of regular extremal solutions of the characteristic 
equation (3.14). So in particular, and in accordance with the notation introduced above, 
the 'BPS substrata' M.\ and A4\ correspond to | and | BPS black holes, all of which 
share the same | super symmetry, respectively (while by .Mi° we denote the stratum of 
general | BPS solutions of vanishing horizon area). In all three cases the £"7(7) invariant 
(){w~^Z) associated to ^ vanishes, and zero, one or three of the eigenvalues of Z AB are 
equal to w in modulus, respectively. These orbits correspond to the various Majorana- 
Weyl spinors of null type which turn out to be slightly more complicated than in the 
case of the null vectors of SO(p, q). For example, the Spin*(12) Majorana-Weyl spinors 
associated to | BPS solutions are pure in the sense of Cartan. 

By virtue of the grading, any linear combination of these charges satisfies ^ 5 = 0, 
and moreover, the commutator of any two of them gives 

= {w n w m -w n w m -2Z nAB Z^ + 2Zf B Z mAB ){a 1 + a 2 ){c?-a l ) G 1 (2) (6.12) 

The generator (a 1 + a 2 ) (a 2 — a\) of 50*(16) leaves invariant the Killing spinor associated 
to the charges and the Killing spinor equation trivially reduces to 

de£ = (6.13) 

We have thus derived the most general | BPS solution of maximal supergravity, which 
depends on 32 independent harmonic functions associated to the independent components 
of the Spin* (12) Majorana-Weyl spinor. Of course in general, some of the ^ n can preserve 
more supersymmetry, in such a way that each black hole can be either a generic | BPS 
black hole carrying one of the 32 linearly independent charges of or a | BPS black 
hole of vanishing horizon area carrying one of the 31 linearly independent charges of Ai\ , 



10 where Ic(SU{2) x Sp{Z) K 14 2 ) is the contracted form of Ic(SU{2) x SU(6)) : 
(SU(2) x Sp{3)) x (((2 <g> 6)+ ffi 14 2 ) (1) © (2 <g> 6)^ © 1 <3) ) 
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or a \ BPS black hole carrying one of the 25 linearly independent charges of M.\, or a 
| BPS black hole carrying one of the 16 linearly independent charges of M\. All these 
solutions can be understood as | BPS multi-black holes of the magic Af = 2 supergravity 
associated to the quaternions [22]. Within the latter truncation, Di is a non compact 
generator of the sl 2 component of sl 2 ©so*(12) C e 7 (_ 5 ), which decomposes e7(-5) as 

e 7( _ 5) l<- 2 > © 32^ © (flfi © so*(12)) (0) © 32« © 1 (2) (6.14) 

exhibiting the fact that these solutions can be understood inside the M = 2 truncated 
theory. 

6.2 ± BPS solutions 

For solutions preserving | supersymmetry the charge state is of the form (using (4.19) 
once again) 

\<jf) = e ^ AB a A a e ^ w + ZAsa A a B + J_ £ ^^ w a A a B a C a 5 ^) |0) (6.15) 

in an appropriate basis, where VLab is a complex-selfdual symplectic form of C 4 and 
Z AB is complex-selfdual with respect with the complementary U(4). The expression in 
parantheses is thus a chiral spinor of Spin* (8); by triality and because of the isomorphism 
Spin*(8) = Spin(2,6) it can be equivalently viewed as a vector of 5*0(2,6). The Killing 
spinors for (6.15) then satisfy 

e£ + e a pQ AB e B = ef = (6.16) 

The associated grading operator of so* (16) is 

Di = ^(Q AB a A a B - Q AB a A a B ) (6.17) 

and we now have 

Di\<g) = 2|#) (6.18) 
The corresponding five-graded decomposition of so* (16) is 

so*(16) = 6 ( - 2) © (4©8)^ 1> © ( [ 1 ©su*(4)©so*(8)) (o) © (4©8)^©6 (2) (6.19) 

One computes that transforms as a Major ana- Weyl Spin* (8) spinor with respect 
with SU*(4) x Spin*(8) = Spin(5, 1) x Spm(2,6), i.e. as a vector of 50(2,6). Indeed 

c 8(8 ) so* (16) = 8^ 2) © (4 © 8_)^ 15 © (6 © 8+) (0) © (4 © 8_)£ © 8$? (6.20) 
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where by 8 and 8± we denote the three inequivalent fundamental representations of 
Spin*(8). In accordance with the grading (6.20) a charge matrix of the form (6.15) 
transforms as 8 + and defines a nilpotent abelian subalgebra R, 2+6 of t$(8)- The corre- 
sponding moduli space is 

R;x 50(2,6) Spin* (16) 

2 50(1,6) 2 (5£/*(4) x Spin(l, 6)) k ((4 x 8)r © 6) 1 J 

The subgroup 5*0(1,6) here is the isotropy group of a timelike vector of 50(2, 6). In 
order for the corresponding multi-black holes solutions to be regular, one must thus 
require the charge matrices in (6.21) to correspond to non-space-like vectors. When the 
norm of the 50(2,6) vector goes to zero, we reach the boundary of M\ corresponding 
to the subspace of \ BPS charge matrices belonging to 

lj9 50(2,6) tJ Spin*(16) 

where we again employ the notation in (6.10). So we conclude that the general ~ BPS 
solutions depend on eight independent harmonic functions associated to the linearly 
independent non-spacelike 50(2, 6) vectors. Note nevertheless that even if all the vectors 
are chosen to be null, the solution only preserves one quarter supersymmetry since, 
although each black hole is | BPS, there are only eight supercharges that are preserved 
by each one of them. As the most general | BPS solutions were possibly understood 
within the magic M = 2 truncated theory associated to the quaternions, the general 
| BPS solutions can be understood within Af = 4 supergravity coupled to six vector 
multiplets. In the latter truncation, Di is the generator of the 50* (8) R-symmetry group 

4 

of the theory, 

so *(8) 6 ( - 2) © (gq © su*(4)) (0) © 6 (2) (6.23) 

such that 

so(8,8)0 (so* (8) ©so*(8)) = 8^' © (6 © 8+) (0) © 8^ (6.24) 
6.3 \ BPS solutions 

For | BPS solutions, the set of 16 supersymmetry charges is associated to the Killing 
spinors 

e l a + e a ^e^ = (6.25) 

where Qij is a (generally complex) symplectic form of C 8 satisfying fijfcfP = Sj (with 
our usual convention Q l i = (Qy)*). In this case the four-dimensional R-symmetry group 



21 



SU (8) permits to rotate fiy to a diagonal basis 

n* = e** jj (6-26) 

but one cannot eliminate the overall phase uj in general. Eq. (4.5) uniquely determines 
the charge matrix as 

\<g) = me-^e^^v |o) (6.27) 

where the overall phase factor e~ 2tbJ ensures that |^) is a Majorana spinor for any choice 
of uj via the identity 

e 2tuJ il[ijQ k i] = e 2 ^ —Eijki m npqQ mn Q P9 (6.28) 

Asymptotically Minkowskian \ BPS black holes thus correspond to uj — 0. The asymp- 
totic central charges are then purely electric modulo an SU(8) rotation. 

The general \ BPS solution is determined by one single harmonic function and the 
given | BPS charge matrix. The relevant grading operator is 

Di = -(%aV - Q ij a iaj ) (6.29) 

and satisfies Di \€) = 4|^). It commutes with the generators of su*(8) and its action 
on the nilpotent generators follows from 

[Di , a? ± tt ij aj ] = ±(a l ± (6.30) 

and thus gives rise to the decomposition of so* (16) 

so*(16) ^ 28 ( ~ 2) © (fl^ © su*(8)) (0> © 28 (2> (6.31) 

such that 

e 8 ( 8 ) e so*(16) ^ 1 ( " 4) © 28 ( " 2) © 70 (0) © 28 (2) © 1 (4) (6.32) 

and = |^ (4) ) G 1 (4) . Let us point out once again that the level ±2 representations 
of SU*(S) are both real, but dual to one another, and correspond to independent real 
combinations of the field strengths, as explained in section 5. 

6.4 Non-BPS solutions 

The results of the preceding section can be summarised by characterising the relevant 
vector space of g © fj*, i.e. g (2) © J) (2) = 32, 8 and 1, for the | BPS, the \ BPS and the \ 
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BPS solutions, respectively. One can now construct non-BPS multi-particle solutions in 
a similar way. Namely, the charge matrices of non-BPS extremal solutions are associated 
to the same graded decomposition (6.32) of e 8 (s) as the maximally supersymmetric ones. 
The only difference is that the associated charge matrices now belong to the grade one 
component of (6.32), which defines an abelian nilpotent subalgebra R 1+27 C es(8)- The 
relevant (a;- dependent) generator H £ so* (16) is thus identical with Di, cf. (6.29). 
The simplest solution of this type is obtained from the state 

|^< 2 >) = (1 + *) (l + -%aV')ie- 2 ^|0) (6.33) 

where the overall phase is determined such that the four-form component of 

(H - 2) |^ (2) ) =1(1 + *) ie~ 2iuJ %Q kl aVaV|0) = (6.34) 

4^-form 8 

vanishes due to the action of the projection operator |(1 + -k). In order to obtain non- 
BPS asymptotically Minkowskian extremal black holes, one must therefore choose u = j. 
The top grade component then becomes associated to a pure | BPS Taub-NUT black 
hole as would be obtained from the level four charge matrix 

=i ne ^i^ |o) (6.35) 

The charge matrix of grade two depends on the associated mass m and a rank two tensor 
Qij satisfying 11 



Qij = fii fc %Q fei n t3 Qij = (6.36) 

such that 

|<T 2 >) = (1 + *) (l + -fiyaV') (m + Q y aW) |0) (6.37) 

In this way one obtains multi-particle solutions depending on 28 harmonic functions. 
Although the symplectic form is not an invariant tensor of SU*(8), the reality condition 

= £lik£ljiZ kl (6.38) 

is nonetheless preserved by SU*(8), and it defines the real 28 representation of SU*(8). 
This is also the representation to which belong the asymptotic central charges 

= —Vtij + Q^ (6.39) 



11 Which means that Qij lies in the 27 representation of Sp(A) C SU(8) which leaves invariant fl 



<.i ■ 
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In contradistinction to the \ BPS asymptotically Minkowski black holes, the non-BPS 
ones are therefore purely dyonic (u = j), up to an SU(8) rotation. 

|<^ (2) ) corresponds to a regular black holes as long as the eigenvalues of Zy are less 
than or equal to m in modulus. The £7(7) quartic invariant §(Z) is then negative, 

0(Z) = ~ + jQ tJ Q l3 + j^Q^QmQ 11 - \(Q ij Q i3 ) 2 + 2Q ij Q^ k Q kl Q H < (6.40) 

In order to spell out the these conditions explicitly, it is useful to find a basis in which 
both Qij and are block diagonal (which always exists thanks to the reality condition 
(6.36)) 





(6.41) 



/ P> \ 

p 2 

p 3 

\ - Pl - P2 - P J 

The three real parameters p t are constrained to lie in the tetrahedron defined by the 
conditions p< < y for i — 1, 2, 3 and — y < p 1 + p 2 + p 3 . In this basis, the £7(7) invariant 
reads 

0(Z) = - (- - p,) (- - p ; ) (— - P ) ( J + P, + P, + p) (6.42) 

and one checks easily that it is negative inside the tetrahedron and that it vanishes on 
its faces. In the absence of NUT charge, one expects the horizon area of such black holes 
to be given by 

Air = 47rv/-0(Z) (6.43) 

Such multi-black hole solutions include non-BPS black holes corresponding to one of the 
28 linearly independent charge matrices in 

K x SU*(8) Sptn*(lQ) 

Sp(4) Sp{4) x 27 1 ' 

(for which the parameters p { lie inside the tetrahedron) as well as | BPS black holes 
corresponding to one of the 27 linearly independent charge matrices in 

SIT (8) c Spm*(16) (6 45) 



(SU(2) x Sp(3)) x (2 ® 6) B R+ x Ic{SU{2) x Sp(3) x 14, 

(for which p t lie on a face of the tetrahedron). Finally, | BPS black holes correspond to 
one of the 22 linearly independent charge matrix in 

SU*(8) c gpin* (16) (6 46) 



(S'C/*(4) x Sp(2)) k (4 <8) 4) R (5pzn(5, 1) x Spm(l, 6)) x ((4 ® 8) R © 6) 
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(with Pi lying on an edge of the tetrahedron) and \ BPS black holes to one of the 13 + 1 
linearly independent charge matrices in 

SU*(8) ^ Spin*(16) . . 

URC — f- fo . K ' (6.47) 



(SU*(Q) x Sp(l)) k (6 <g> 2) R 5C/*(8) k 28 

(corresponding to p ; lying on a vertex of the tetrahedron). This set includes the charge 
matrix (6.35) lying in the component of grade four 1 (4) . 

Note that these solutions are completely different from the BPS ones. To see this, 
we note that the decomposition of e§(8) associated to a \ BPS charge matrix within the 
graded decomposition associated to either the | BPS solutions or the non-BPS extremal 
solutions lead to equivalent decomposition in terms of representations of SU (2) x S77*(6). 
The charge matrix in the neighbourhood of this charge which define ~ BPS black holes 
which do not interact with the | BPS black hole associated to the latter, can be obtained 
by acting with either so* (12) or su*(8), for | BPS and non-BPS solutions, respectively. 
The action of so* (12) generates new components in the 15 of 577* (6) whereas the action 
of su*(8) generates new components in the (2 ® 6)^ of SU(2) x S77*(6), which shows 
that they are inequivalent. 

Acknowledgments: We are grateful to Boris Pioline and Kelly Stelle for discussions 
and comments. 
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